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Abstract 

Boundary-induced transport in particle systems with anomalous diffusion 
exhibits rectification, negative resistance, and hysteresis phenomena depend- 
ing on the way the drive acts on the boundary. The solvable case of a ID 
system characterised by a power-law diffusion coefficient and coupled to two 
particles reservoirs at different chemical potential is examined. In particular, 
it is shown that a microscopic realisation of such a diffusion model is provided 
by a 3D driven lattice- gas with kinetic constraints, in which energy barriers 
are absent and the local microscopic reversibility holds. 
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Introduction. - Macroscopic motion generally results from the action of nonzero macro- 
scopic forces. Ratchets are systems which are able to develop a directed motion in the ab- 
sence of macroscopic forces. Since the early days of kinetic theory, devices of this sort have 
been used as a means to probe the statistical nature of the second law [I]]. A celebrated 
example is the Smoluchowski-Feynman ratchet and pawl machine [|2],|3|]. In the last decade, 
the interest in ratchet systems has been mainly motivated by motor proteins and new sep- 
aration techniques [Q. According the Curie principle, a directed transport may arise even 
in absence of a macroscopic bias provided that both parity and time reversal symmetry are 
broken. Closely related is the non-monotonic behaviour of the particle current in response 
to a driving force, known as negative incremental resistance (NR). Non-linear transport 
properties are a crucial ingredient in the complex behaviour of many biological and artificial 
systems. Typical NR rectyfing devices are the tunnel diode and the sodium channel. There 
is interest in understanding the microscopic origin of such a behaviour (which in tunnel 
diode is quantum mechanical) and to provide stochastic analog |5||J . Cecchi and Magnasco 
have suggested a purely geometric mechanism in which the time scales involved in the par- 
ticle motion do not follow the Arrhenius- Kramers law but rather depend on the existence of 
purely "entropic" barriers 0. 

In this note we show that similar non-linear transport phenomena may occur in a kinetic 
lattice-gas model in which there are no energetic barriers and no local breaking of the 
detailed balance (time- reversal symmetry). We first consider transport properties in a ID 
solvable model of boundary- driven system with a power-law vanishing diffusion coefficient. 
The anomalous diffusion coefficient induces a non-linear relation between the particle current 
and boundary densities, which in turn is responsible for rectification and negative resistance. 
For the particular case we consider here these features can be analytically investigated. We 
then show that a microscopic realization of such restricted diffusion model is provided by a 
3D boundary-driven lattice gas with reversible kinetic constraints which is coupled to two 
particle reservoirs. 
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The diffusion model. Consider a transport process in a slab of size 2L which is 
described by the one-dimensional diffusion equation: 



dp d 
dt dz 



(i) 



where p(z,t) is the local particle density, and \z\ < L. The diffusion coefficient, D$(p), 
vanishes at a critical threshold density, p c , as a power-law with an exponent (p > 0: 



D^ P ) = D o (i + 0)(p c - P y 



(2) 



The system boundaries, located at z — ±L, are in diffusive contact with two particle reser- 
voirs at chemical potential /i± which keep the boundary densities at 



p{±L,t) =p±, Vt > 0. 



(3) 



When p + = p_ < p c , the characteristic relaxation time is finite and the system attains an 
equilibrium state characterised by a flat profile. When p± = p c , the system approaches the 
critical density by a power-law and the breaking of time-translation invariance ensues j9]. 
The model was indeed introduced in ref. with the purpose of understanding aging in 
a kinetically constrained lattice-gas Here we consider the non-equilibrium stationary 

properties of a system with p + ^ p_ and p + , p_ < p c . In this case the relaxation time is 
still finite for any finite L, and the steady state density profile, p(z), is easily computed 



p c — p(z) = (L a + — z a J) 



1 + 



(4) 



where the constants a± are determined by the boundary condition 

I 



a± = — 
± 2L 



(Pc-P-) 1+ ^±(Pc-P + ) 1+0 



(5) 



The particle current is then obtained as J = D(p)d z p, which gives 



J(p+,p- 



Do 
2L 



(pc-p-) 1+ Mpc-p + ) 1+0 



(6) 



The expression of the current has two interesting features: it is non-linear and does not 
depend only on the single variable Ap = p + — p_. In the limit p± p c the density profile 



is linear, and the Fick's law J ~ Ap is recovered in agreement with the linear- response 
theory. At high enough density, however, this is not the case and more interesting transport 
phenomena emerge. In the following we explore the implication of Eq. (|6|) in the non-linear 
regime for some relevant cases. 

Rectification. To begin with we consider, as in JTIJ, a boundary potential Ap(t) = 
p+(t) — P-{t) which is a periodic asymmetric step function of time with zero average over 
the period r: 

1 - to/t , t G [0, r ] ; 
Ap(t) = \ (7) 

-t /t, t e [t q ,t] . 

The average current over the period is: 

Jav = - r J\p + (t),p-(t)}dt. (8) 
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In Fig. m we show a plot of J av vs the asymmetry parameter a — t/tq — 1, for the case in 
which p-(t) = for t G [0, r ], and p+(t) = for t G [r ,r]. For a = 1 the potential is 
symmetric, Ap(ro + £) = — Ap(t), and no net average current can flow through the system; 
the rectification effect, i.e. a non-zero average net current, occurs for any a^l. As in other 
ratchet systems, there is a certain value of a for which optimal pumping condition exists. 
The current direction is determined only by a: for a > 1 the current is negative and such 
that: J av (o;) = — J av (l/a). However if the density of one boundary, say p_ is kept fixed 
to a non-zero value, there is no current inversion: J av is always negative (positive) regardless 
the value of the asymmetry, a, see Fig. [TJb. In this case, the net current is non-zero even for 
a symmetric potential Ap(t), provided only the diffusion is not normal, ^ 0. This is not 
in contradiction with the Curie principle, as parity is explicitly broken: p(—L,t) ^ p(L,t), 
while the boundary dissipation breaks the time reversal symmetry. More generally, one may 
also consider a symmetric potential that changes adiabatically according to suitable paths 
in the space of variables (p + ,p_), and finds similar rectification properties. 
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Negative resistance and hysteresis. - To further explore the nonlinear transport regime 
we looked for situations in which an increasing driving force leads to a decreasing particle 
current. To keep things simple, we consider the case in which the ratio of the boundary 
densities, 5 = p~/p + , is fixed, and study the stationary point of the Eq. (|6|) at constant 
5. This is done by the method of Lagrange multipliers, which gives the maximum of the 
current 

/ = f£ \1 J. (q) 

at 

p max l _ 5 ± max 1 _ 

— = JM , ^— = 5 tw- (10) 

Pc 1 _ § 4- Pc 1 — (W 

Increasing the driving force above this value leads to a decreasing particle current. The non- 
monotonic behaviour of the current is shown in Fig. ^| as a function of the drive Ap = p + —p- 
for different 5. In the limit of small density gradient the Fick's law is correctly recovered 
(the dotted line in Fig. The higher the 5 the smaller the range in which the linear 
relation holds. The negative resistance found here is a consequence of the fact that the 
current depends in a non-linear way on both reservoir densities, which in our model follows 
specifically from a power-law vanishing diffusion coefficient. This last feature is commonly 
observed, e.g. in colloids and hard-sphere systems near their random close-packing limit. In 
the next section we shall see how, at microscopic level, negative resistance and rectification 
may occur in the absence of any energetic barriers to the particle motion, through a purely 
"entropic" mechanism as has been suggested in ref. |J. These phenomena are intimately 
associated with the possibility of having a hysteretic response. If we consider a loop in the 
plane (p + , p_) it can be shown that the particle current responds to the driving force by 
following a clockwise hysteresis cycle (see Fig. |3|). Finally, we mention that similar non- 
linear transport properties also appear in the presence of a divergent diffusion coefficient, 
< 0. The related diffusion equation was studied by Carlson and coworkers in an attempt 
to provide a continuum description of self-organising critical systems JIT]. 



The driven lattice glass. We now show that a microscopic realization of the above 
diffusion model does indeed exist. It is provided by a kinetically constrained lattice-gas, 
boundary driven in a non-equilibrium stationary state by a chemical potential difference. 
The bulk dynamic of the model (defined on a cubic lattice) consists in moving a particle to a 
neighbouring empty site if the particle has less than 4 nearest neighbouring particles before 
and after it has moved, consistently with detailed balance @. The boundary dynamic mimics 
the contact of the system with a particle reservoir at chemical potential p±. It is simulated 
according the usual Montecarlo rule: if a randomly chosen site on the layer is empty, a new 
particle is added; otherwise, if the site is occupied, the particle is removed with probability 
e~ M± (p± > 0, we set k-^T = 1). The global effect of the reservoirs is to fix the boundary 
densities at two different values p + and p_ , driving a current through the system. The aging 
dynamics of the system when the two reservoirs are at same chemical potential was first 
investigated in |8|||. As the density gets closer to the threshold value p c ~ 0.88, the diffusion 
coefficient approaches zero as a power-law, Eq. (|2|), with = 3.1 0. For our purposes it 
is sufficient to show that close to the threshold p c and in the presence of two reservoirs, 
the system approaches a steady state characterised by a density profile which is exactly 
described by the Eq. ([|). In Fig. [|the density profile obtained in a Monte Carlo simulation 
is compared with the one predicted by the anomalous diffusion equation using the above 
values of p c and <ft. There is excellent agreement between the two. The small discrepancy 
observable near the higher density edge is a finite-size effect which tends to disappear as the 
thermodynamic limit L — > oo is approached. In Fig. f| we show for comparison the numerical 
density profiles of both the usual boundary driven 3D lattice-gas (obtained by removing the 
kinetic constraints), and that predicted by the normal diffusion equation. These results 
suggest that the non-linear nature of the density profile, and consequently the non-linear 
transport, is essentially determined by the presence of blocked configurations induced by the 
kinetic constraints. 
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Summary and conclusion. - In this paper we have shown that an analytically solvable 
model of boundary-induced transport exhibits rectification, hysteresis and negative resis- 
tance phenomena. We have also shown that such features may occur in driven lattice-gas 
models where local detailed balance holds and dissipation is only forced on the boundary In 
particular, we have presented numerical results supporting the hypothesis that the diffusion 
model represents the hydrodynamic limit of a driven lattice-gas with kinetic constraints. 
The efficiency of the system in presence of an external load and the interplay of the bound- 
ary driving force with a spatially varying potential may also be of interest. Since there are 
no energetic barriers to the motion of particles this model provides another example of a 
ratchet system based on the mechanism of entropic trapping ||. The relationship between 
the macroscopic diffusion model and the microscopic lattice-gas sheds some light on the 
non-linear nature of transport properties. On one hand, they are due to the non-linear de- 
pendence of the density profile on both boundary densities, not simply on their difference. 
On the other hand, the non-linear density profile in the driven lattice-gas follows from the 
presence of blocked configurations. This suggests that the steady state transport properties, 
like some features of aging systems, are essentially determined by an extensive entropy of 
blocked configurations. It is tempting to speculate that even in such cases the invariant 
dynamic measure could be defined in terms of a suitable generalisation of the Edwards 



hypothesis 12 
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FIG. 1. Rectification. Average current, J av , plotted against the asymmetry parameter 
a = t/to — 1 for a boundary periodic potential with zero average over the period [0, r] (in both 
figures we set r = 1). a) p~(t) = for t € [0,to], and p+(t) = for t € [to,t]. b) The curves with 
J av > (J av < 0) correspond to a fixed p + (p~). In both cases there is a value of a for which 
optimal pumping condition exists. 
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FIG. 2. Negative resistance. The particle current J(p+,p-) is plotted versus the potential 
p + — p- for different value of the ratio 5 = P-/ p+ = 0.3 , 0.4 , 0.5 , 0.6 (from top to bottom). The 
dotted line represents the Fick's law, which is recovered in the limit of small density gradient. 
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FIG. 3. Hysteresis. Response of the current, J, to a cyclic change of the potential 
Ap = p + — p_. The loop in the plane (p+, p_) is carried out in three steps: starting from 
p + = P- = p a the value of p + is increased up to pb keeping constant p_ = p a (upper curve); then 
p_ is increased from p a to pb while p + is fixed to pb (lower curve); finally p + and p_ are decreased 
to p a keeping p + = p_. We set p a = and pb = 1. The continuous line corresponds to the case 
p c = 1 and cf) = 3.1 (anomalous diffusion); while the dashed straight line - absence of hysteresis - 
to normal diffusion (0 = 0). 
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FIG. 4. Density profile in a boundary-driven lattice glass coupled to two particle reservoir 
at density p + = p(L) = 0.85 and p_ = p(—L) = 0.75 (squares), with L = 160 and transverse 
surface 20 2 . The continuous smooth line represents the analytic profile predicted by the diffusion 
equation, Eq. (Q). Also shown, for comparison, are the linear profile corresponding to a normal 
diffusion coefficient (dashed line), and the numerical simulation data obtained by removing the 
kinetic constraints (stars). 
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